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The behavior of a heavy tagged intruder immersed in a bath of particles evolving under ballistic 
annihilation dynamics is investigated. The Fokker-Planck equation for this system is derived and 
the peculiarities of the corresponding diffusive behavior are worked out. In the long time limit, the 
, intruder velocity distribution function approaches a Gaussian form, but with a different temperature 

from its bath counterpart. As a consequence of the continuous decay of particles in the bath, the 
mean squared displacement increases exponentially in the collision per particle time scale. Analytical 
results are finally successfully tested against Monte Carlo numerical simulations. 

^sO ■ PACS numbers: 51.10.+y,05.20.Dd,82.20.Nk 

Q ■ I. INTRODUCTION 

<d ; 

In recent years, there has been some interest for systems where particles annihilate ballistically [H, 0, H, 0, S S 0, II] ■ 
In these studies, the model considered consists of an ensemble of hard particles which evolve freely until a binary 
encounter, which leads either to the annihilation of the colliding partners with probability p, or to an elastic collision 
with probability 1 — p. For this probabilistic annihilation model, there are no collisional invariants, and numerical 
simulations have shown that for a broad class of initial conditions, the system reaches an homogeneous state in which 
all the time dependence of the one particle distribution function is encoded in the density and temperature (defined as 
the second velocity moment of the distribution function) 0, Q . This is the so-called "Homogeneous Decay State" . Such 
a behavior resembles the one of granular fluids (see Q and references therein) where, if the system is stable, it evolves 
into an homogeneous cooling state, in which all the time dependence is borne by the granular temperature (in this case 
the density is conserved) Il0|. For the annihilation model, the hydrodynamic equations have been derived using the 
Chapmann-Enskog method [ll| by the usual assumption of the existence of a "normal solution" , whose space and time 
dependence occurs only through the hydrodynamic fields [f| . Recently, the hydrodynamic equations linearized around 
I \ the homogeneous decay state have been derived relaxing such an assumption Nevertheless, it must be assumed 
that there is scale separation, i.e that the spectrum of the linearized Boltzmann collision operator is such that the 
■^j- ' eigenvalues associated to the hydrodynamic excitations are separated from the faster "kinetic eigenvalues" . Although 
. this property is valid for elastic collisions fla|. it has not been proven for the probabilistic ballistic annihilation model 
^^O ' in general, but only for Maxwell molecules and for p smaller than a given threshold [12| . 

The objective in this paper is to study the simplest transport process in this system in which we can rigorously 
y—i ■ prove that there is scale separation. We will consider a tagged particle in a fluid in the homogeneous decay state, 
but collisions between the tagged particle and the particles of the fluid will be always elastic. The equation for the 
tagged particle is the Boltzmann-Lorentz equation (l3l . [l5| which depends on the one particle distribution function 
of the bath. In the limit of asymptotically large relative mass for the tagged particle, this equation reduces to a 
Fokker-Planck equation which depends on the time-dependent density and temperature of the bath. Due to the 
structure of this equation, we can prove that there is scale separation and that, in the long-time limit, the velocity 
distribution function of the tagged particle approaches a Gaussian distribution but with a temperature that differs 
from that of the bath. A similar breakdown of equipartition has been reported for a heavy particle in a granular bath 
fl7l. HH. [l9L I20I I2l1] . a problem that can be mapped onto an elastic situation b 22], at variance with the situation under 
scrutiny here. We also study the diffusion of the heavy particle and identify the diffusion coefficient as a Green-Kubo 
formula in terms of the velocity autocorrelation function. Finally, we perform Monte Carlo numerical simulations in 
order to test our theoretical results. 
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II. FOKKER-PLANCK EQUATION 



We consider a tagged particle of mass m and diameter a immersed in a low-density gas. This gas is composed of 
hard spheres or disks of mass m g and diameter a g which move ballistically until one particle meets another one; such 



2 



binary encounters lead to the annihilation of the colliding partners with probability p or to an elastic collision with 
probability 1—p [2, [2, [H, 0, H, 3 ■ Collisions between the particles of the gas and the tagged particle are always elastic. 

A. From Boltzmann-Lorentz to Fokker-Planck 

The evolution equation for the probability density F(r, v, t) of the tagged particle is the Boltzmann-Lorentz equation 

PHI 

^+v-VjF(r,v,t) = J[r,v,t|F,/], (1) 

where the collision operator is given by 

J[v,v,t\F,f] = a*- 1 fdv 1 fd&Q( g -&)(g-&){F(T,v*,t)f(v,y* 1 ,t) 

- F(r,v,t)/(r, Vl) t)}. (2) 

Here /(r, v, t) is the distribution function of the particles in the gas, d is the space dimension, g = v — vi is the 
relative velocity, O is the Heaviside step function, it is a unit vector pointing from the centre of the gas particle to 
the centre of the tagged particle at contact, and ctq = a+ 2 9 ■ The precollisional velocities v* and are given by 

2A 

v * =v- — ^(g-o-)*, (3) 
Vl =v 1 + ^-^( g -<r)<T, (4) 

with A = rrig/m the (gas/tagged particle) mass ratio. 

We shall consider that the gas is in the homogeneous decay state, so its distribution function has the scaling form 

1 

/ff(vi ' t) = <w Xff(ci) ' Cl = wr () 

/ 2T ft) \ 

where n g (t) is the number density of the gas, v g (t) — t ^ J is the thermal velocity of the particles in the gas 

and xh is an isotropic function depending only on the modulus c = |c| of the rescaled velocity. It can be seen that 
the homogeneous density and temperature obey the following equations 0] 

" -pu g (tX n n g (t), (6) 



dt 

■A 

-pv g (t)C T T g (t), (7) 



ar g (t) 



dt 

where we have introduced the collision frequency of the corresponding hard sphere fluid in equilibrium (with same 
temperature and density) 

d-l 



i 1 / 2 (d + 2)r(d/2)' 



».®= gW I, ? 9 TT^TTT^ - (8) 



Here the dimensionless decay rates £n and Cr are functionals of the distribution function and are approximately 
known in the first Sonine approximation 0, Q , see Appendix [XJ Equations (J6j) and |[7J can be integrated to obtain 
the following power laws for the decay of the density and temperature 

n g (t) = n a (Q)[l + v g {0)p(C n + C T /2)t]-*&te , (9) 
T g (t) = T 9 (0)[l + ^(0MC„ + CT/2)r^fe. (10) 

As a consequence, we get n g T g ^ 2 oc t~ l , a simplified form of a scaling relation common to all ballistically controlled 
processes fig . 
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We next study the evolution equation for the tagged particle in the limit of large relative mass for the tagged 
particle. In the limit A <C 1, it is possible to expand the collision operator J[r, v, t\F, /] in powers of A. In Appendix 
lElit is shown that the leading order is 

J[r,v,t\FJ] ~ A • [A(v)F(r,v,t)] + : [N(v)F(r,v,t)}, (11) 

where 

A(v,t)= 7 (t)v, N(v,t) = Vy(t)I, (12) 
and 7 is the second order unit tensor. The definitions of 7 and 7 are respectively 

7 (f) = le [n g (t),T g (t)]a(p), (13) 

7(f) = le [n g (t),T g (t)}a(p)Kp)^-, (14) 

The friction coefficient 7e (i) is the same as for elastic bodies, and appears here as a function of the time-dependent 
density n g (t) and temperature T g (t) 



d-l 



with a(p) and 6(p) functionals of the distribution function of the bath which depend only on the parameter p 

a{p) = r((d+V)/2) / rfc i^(ci) Cl , (16) 

KP) = —1 Jdc lXH ( Cl W (17) 

In Appendix [Bl it is shown that the two terms on the right hand side of equation (jlip are both of order n g v g cr ~~ A, 
while the other contributions in the Kramers-Moyal expansion are at least of order rigVgCTQ -1 A 3 / 2 . In the same 
Appendix, the expressions for a(p) and b(p) are evaluated to first order in a Sonine expansion. 

Taking into account the approximate expression for the collision operator, Eq. (|11|) . it is possible to write the 
Boltzmann-Lorentz equation as a Fokker-Planck equation for asymptotically small A 



d 



F{v,w,i)= le {i)a{p)^ 



Hp) 



Tgit) d 

m (9v 



F(v,v,t). (18) 



As in the inelastic case, the Einstein relation is violated due to the fact that the distribution function of the bath 
is not Maxwellian [3 [23|, [24|, H^], which in turn implies that b(p) ^ 1. On the other hand, if we suppose that the 
velocity of the tagged particle obeys a Markov process and write the corresponding Fokker-Planck equation, in terms 
of the jump moments, limAt— yQ ( Au) /At and limAt^o(Av 2 )/A£, we obtain exactly Eq. I|18p. Here (...) means average 
over different noise (bath) realizations. 

B. Coarse grained fields and relevant scales 

We now focus on the study of the hydrodynamic fields of the tagged particle with the aid of the Fokker-Planck 
equation, Eq. (|18p . We define the mean velocity and the temperature of the Brownian particle as 



u(t) = JdrJdwF(r 7 v 7 t), (19) 

~T(i) = /dr/dv~m(v-u) 2 F(r,v,t). (20) 

Taking moments in the Fokker-Planck equation, we obtain (see Appendix |C|) 

du(t) 



dt 
dT(t) 



7e(i)a(p)u(t), (21) 
2 le {t)a{p)[T{t)-b{p)T g {t)}. (22) 
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As the function j e is a known functional of the gas density n g and temperature T g , equations pTj) and (f22|) can be 
integrated, which yields 



and 



u(i) = u(0) [1 + i/ fl (0)p(C n + Cr/2)i] <2C -+^ 



(23) 



1 - 6 

T(Q)- 



&(p)r a (o) 



1 -e 



[1 + i/ fl (0)p(Cn + Cr/2)t]~^R53 



In the above equation, we have introduced the dimensionless coefficient e 

vCrv g {t) _ V2d( T 



2a{p) le (t) 2(d + 2)a(p) 



<?0 



d-1 



A' 



(24) 



(25) 



that is not necessarily a small quantity. 

As can be seen in Eq. (|24p . the behavior of the temperature depends strongly on the value of e. If e < 1 the first 
term of equation (|24p dominates in the long time limit and the temperature of the tagged particle asymptotically 
decays with the same power as the temperature of the gas (see equation (fTU)) '). As a consequence of (|2"4")l we have 



T(t) b(p) 



lim 

t^ooT g (t) 1-e 



e < 1. 



(26) 



On the other hand, if e > 1 the second term of equation (|24[) dominates in the long time limit and the temperature 
decays slower than the gas temperature. One can understand this behavior as follows; the parameter e is essentially 
the quotient between the cooling rate of the gas and the relaxation rate of the tagged particle's temperature. If the 
former is smaller than the latter, the tagged particle's temperature is eventually slaved by T g due to the second term 
of the Eq. (|22|) . In the reversed case, the tagged particle's temperature evolves independently in the long time limit 
with a cooling rate slower than that of the gas. It should be emphasized here that in the expansion made in the 
previous section, we implicitly assumed that T/T g remains finite, because the coefficients A and N were expanded 
in powers of A 1 / 2 ^ (see Appendix [Bl) . Hence, the Fokker-Planck equation for e > 1 is restricted to a time window 
in which t2- is small enough. In the following, we will only consider the case in which the Fokker-Planck equation is 
valid for all times, i.e. the double limit 



A^O, p^O, 



e<l, 



e cx 



do 



P_ 
A' 



(27) 



where the requirement of small p stems from A < 1 and e bounded from above. Hence, in order to be consistent with 
this limit, we can substitute in the Fokker-Planck equation, Eq. ([15]). the values of the coefficients, a(p) and b(p) by 
their elastic limits 



lim a(p) = 1, 

p— >o 



lim b(p) = 1, 



so that 



d 



F(r,v,i) 



Tg(t) 8 



<9v 



F(r,v,t). 



(28) 



(29) 



This equation is formally identical to the one obtained for an elastic gas [l5j], except for the fact that the density 
and temperature of the gas depend on time. This is a consequence of the elastic limit to which we are restricted. In 
general, the coefficients a(p) and b(p), Eqs. (fT^|) and (fT7)) . differ from unity due to the non Maxwellian character of 
the distribution function of the bath, and Einstein relation is violated. 

In order to analyze this equation, it is convenient to introduce the dimensionless time scale, t* , proportional to the 
number of collisions of the tagged particle 



t* = (l-e ) / dt' le {t'), 
lo 



(30) 
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where cq is defined by substituting (t(p) / 'a>(p) m (EH) by its p — > limit: 



The dimensionless time scale t* is related to the real time t by 

t* = 1^—^. J_ T log[l + i/ fl (0)p(C„ + Cr/2)t]. (32) 
In this time scale, the evolution of the mean velocity and temperature of the tagged particle are particularly simple 

u(t*) = u(Q)e'^, (33) 



and 



T(f) _ T(Q) „ t . , 1 . „ t . 



e~ 2t + - (l-e~ M ). (34) 



T g {t*) T g (0) l-e, 

Such predictions will be compared against numerical simulations in section IIVI Let us also introduce the scaled 
distribution 

1 



where 



"(r,v,t) = W) F-(r,V,n V = — , (35) 



The function v e (t) is introduced because with these definitions we have 

"2T(t) ll/2 



Ve(t) 



m 

in the long time limit. In these variables the Fokker-Planck equation (|29|) reduces to 



(37) 



JL + £ (t*)v* ■ V ) F*(r,v*,0 - A FP (v*)F*(r,v*,t*), 



where we have introduced the standard homogeneous Fokker-Planck operator 



<9v* \ 2 <9v 
and the function proportional to the mean free path 

v e (t) - ■ - ' 1 



to(f) = , . . = dT{d/2) s i A' 1/2 f ^ Kfi*)^- 1 ]- 1 . (40) 

V (l-eo)7e(t) (1-60)3/24^^ V^oy 1 L 91 ^ 9 J K ' 

Taking into account the definition of e, Eq. (125|) . we can write explicitly £o as a function of the t* variable as 

e (t*)=to(py' tm , e*= 2Cn€0 (41) 

Ct(1 - eo) 

To sum up, we have obtained the evolution equation for the distribution function of a tagged particle in a bath of 
particles which annihilate, in the limit where the tagged particle is much heavier than the particles of the bath. There 
are some points in common with the elastic case, but also some important differences. The homogeneous operator, 
whose spectral properties are well-known [ID, is exactly the same, but the flux term is weighted by a function 
depending on time and that diverges in the long time limit. This will have important consequences in the study of 
diffusion as we will see in the following section. Moreover, the equation is not valid for all values of the probability p 
of annihilation in the bath and masses of the tagged particle but, as already mentioned, is limited to the double limit 
of Eq. l[2"7j). in which e < 1. 
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III. LONG-TIME LIMIT SOLUTION OF THE FOKKER-PLANCK EQUATION 

In this section we investigate the long time behavior of the solution of the Fokker-Planck equation, Eq. (|38[> . 
starting with an arbitrary initial condition. The objective is to study if the tagged particle reaches some scaling state 
in the long time limit and also to analyze how the particle diffuses. 

A. Evolution towards a scaling form 

As the Fokker-Planck equation is linear, it is convenient to work in the Fourier space. The Fourier component of 
the tagged particle distribution function is defined as 

F k (v*,t*) = Jdre- lkr F*(r,v*,t*), (42) 

so that Eq. ([38]) yields F k (v*,i*) 

£pFk(v*, O = [A FP (vl ~ i*o(t*)k ■ v*]F k (v*,i*). (43) 

The spectrum of the operator App(v*) — i£ (t*)k ■ v* is known [l3|, [l5j]. The eigenvalues are 

d 

MM) = --[M (*)] 3 -$>,-, (44) 

where we have introduced the vector label n = (rii, . . . , n>d), with possible coordinate values rii = 0, 1, 2, . . . oo. 
Hence, for any initial condition, all the k-Fouricr components decay and only the k = remains. Moreover, as the 
eigenfunction associated to the vanishing eigenvalue is the Maxwellian distribution [l3|, [lj| 

= ^- v '\ (45) 

we obtain 

F(r,v,t) ^ ^ X m(v*), (46) 

in the long time limit. 

As a consequence, the tagged particle distribution function approaches a scaling form similar to ([5]) for the gas, but 
with a different temperature (see (f2l)|)). In this regime the cooling rates of the bath and of the tagged particle are 
the same and the temperatures are proportional. The situation is similar to that of an clastic particle in a bath of 
inelastic grains [HI, E^- Nevertheless, there is an important difference: in the inelastic case it has been proved that 
there exists an exact mapping with an elastic system. On the other hand, in our problem such a mapping fails due 
to the flux term which explicitly depends on time. 

B. Characteristics of diffusive motion 

Our objective is to study the evolution equation for the density of tagged particles, n(r,t) — J dvF(r,v,t) in a 
"macroscopic" scale, i.e. in a long time and length scale compared to the microscopic ones. The microscopic time 
scale is defined by the slowest kinetic modes of App, i.e. the modes with a single non vanishing component, labeled 
by rii = &ij for a given value of j in [1, d]. The microscopic length scale is defined by the mean free path of the tagged 
particle which is proportional to £o(t*)- The starting point will be the Fokker-Planck equation for F k , Eq. (|4"3"|) . As the 
generator of the dynamics, the operator App(v*) — i£o(^*)k • v*, does not commute with its time derivative, it is not 
possible to write the general solution of equation (|43[) in terms of the initial condition in a simple way. Nevertheless, 
it is shown in Appendix [Dl that in the hydrodynamic limit, i.e. k£o(t*) <C f for all the time evolution and t* f, a 
closed equation for the Fourier component of the density, n k = J dv*i ; k (v*, £*), is obtained 

d -^l = -D Q [kl Q (t*)]W{t*), (47) 
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where the diffusion coefficient is 

*> = 2(lW (48) 

This asymptotic behavior can be evaluated, taking advantage of the scale separation (i.e. the mode with n = is 
isolated from the other modes). From equation (|47[) we can derive the evolution equation for the density 

Mr,**) _ n ,2^*72 



dt* 



= D £' (t*)V 2 n(r,t*). (49) 



This is the equation we were looking for and it is only valid in the "macroscopic" time and length scale. If we 
transform this equation to real time with the aid of formula (|30p we obtain 



= D (t) v 2 n(r, t), D(t) = C ^ — — — £> e (i), (50) 

dt (1 - c )[Ct + (2C« - Cr)eoJ 

where D e (t) = 2v 2 (t)/^ e {t) is the same as the diffusion coefficient for elastic collisions except that it appears here as 
a function of the time-dependent gas temperature and density. As can be seen, for our system the diffusion coefficient 
D{t) is far from being a trivial generalization of the elastic diffusion coefficient. 

Now let us focus on the predictions of our diffusion equation. To this end, we introduce the mean square displacement 

(r 2 (t*)) = I drr 2 n{v,t*). (51) 



If we consider an infinite system, we obtain from equation (|49|) 

JL(r 2 (t*))=2dD £ (t*) 2 , (52) 

that will only be valid in the long time limit. It is straightforward to integrate equation (|52ll . taking into account the 
explicit formula for £o(t*), Eq. (|4"Tj) . This gives 

(r 2 (O)=^A4(0) e2£ l~\ (53) 

or in real time 

(r 2 (t)) = {[I + ^(0)KC« + Ct/2)*]^ l} . (54) 

As can be seen in Eq. (f53|) . the diffusive behavior is completely different from its elastic or even inelastic counterparts, 
where it was found that the mean square displacement is proportional to the number of collision per particle [la . . 
The fact that the bath is loosing particles significantly affects this dynamics, and the mean square displacement 
increases exponentially in the collision per particle scale. As we will see in the following section, simulation results 
agree well with our theoretical prediction. Roughly speaking, the exponent 4£ n /(2£ r! + (r) is close to id/ {Ad + 1) 
(approximately 1.77 in two dimensions, and 1.84 in three dimensions). 



C. Diffusive behavior : an alternative derivation 



In the remainder of this section, we show that, under plausible hypothesis, it is possible to re-derive "a la Einstein" 
the formula for the mean square displacement, Eq. (|53|) . This derivation has the merit of leading to a Green-Kubo 
like expression for the diffusion coefficient. We start by writing the mean-squared displacement as 

(r 2 (t))= fdt' f dt"{V{t')-\{t")). (55) 
Jo Jo 

Here, the position and the velocity of the tagged particle, r(t) and V(i), are considered as a stochastic process and 
(. . . ) denotes an ensemble average over different trajectories. Let us change the variables from t — > t* and let us also 
introduce the scaled velocity 

w <**> 55 WY (56) 
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Figure 1: (Color online) Evolution of the temperature ratio as a function of the number of collisions per particle r for a system 
with p — 0.01 and two values of the tagged particle's mass, m = 20m fl (A = 1/20) and m = 200m s (A = 1/200). The dashed 
line is the theoretical prediction given by Eq. (|34|l . 



where v e (t) is defined in 



With these definitions we have 



(r 2 (t*)) 



1 



(1 



.., , d Sll -\ Sl ) [ ds 2l - 1 (s 2 )(V(s 1 )-V(s 2 )) 
eo) 2 Jo Jo 



= *§(0) f dai f d S2 e £ * (si+S2) (w( Sl )-w( S2 )), 
Jo Jo 



(57) 



where we have used the definitions of t* and £o(0), Eq. (f50"| and Eq. (|4T|) . Now, if we assume that the tagged particle 
is in the scaled regime, i.e. the temperature is T g (t*)/(1 — eo) and that the correlation function (w(si) ■ w(s 2 )) is 
a function of si — s 2 , by integrating in the new variables S = (si + s 2 )/2 and s = Si — S2, the following relation is 
obtained 



(r 2 (r)> 



<(0) 



ds(w(s) • w(0))( 



(58) 



in the long time limit. This formula is the generalization of the Einstein formula for the diffusion coefficient of a heavy 
particle in a fluid in the homogeneous decay state. It relates the asymptotic behavior of the mean-square displacement 
with the time integral of the autocorrelation function of the velocity weighted by the exponential e _e s . So far we 
have considered the Fokker-Planck equation as the equation for the one-time probability distribution function. If 
we assume that the velocity of the tagged particle is a Markov process, then the Fokker-Planck equation is also the 
equation for the conditional probability and we can evaluate easily the correlation function (w(si) • w(s 2 )). Taking 
into account Eqs. (|33|) and (f36j) . we have 



(w(si) • w(s 2 )) 



-|«l- s 2 



(59) 



By substituting this formula into Eq. 
written in the Green-Kubo form 



we re-derive Eq. (|53|) with the same diffusion coefficient Do, that can be 



Do = -. ds{w(s) ■ w(0))e- 
dJo 



(60) 



IV. DIRECT SIMULATION MONTE CARLO RESULTS 



The objective of this section is to put to the test the main results of the previous sections by means of the direct 
simulation Monte Carlo method (DSMC). More precisely, we will analyze the temperature and mean velocity evolution, 
together with the tagged particle diffusion. We have performed DSMC simulations of a system of N g hard disks of 
mass m g and diameter a g which annihilate with probability p or collide elastically with probability 1 — p evcrytimc 
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Figure 2: (Color online) Scaling function $ defined in the main text as a function of reduced time t* for different systems with 
p = 0.1 and p — 0.01 and different values of the taeeed Darticle's mass. The dashed line is the theoretical prediction of Eq. 

El. 
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Figure 3: (Color online) Stationary values of the temperatures ratio for a system with p = 0.01 and different values of the 
tagged particle's mass A = m g /m. Symbols are for the Monte Carlo data and the dashed line shows the long time limit of Eq. 



two particles meet. Bird's algorithm [26j has been used. The parameters in all the simulations were m g 



1, Co 



1, 



N g (0) — 10 and T g (0) — 1. We have considered only one tagged particle in each simulation, that collides elastically 
with the surrounding bath. The diameter of this particle has been set to unity (a = 1) and we have varied the value 
of its mass m. The values of the probability of annihilation p of the particles in the bath have been p = 0.1 and 
p = 0.01, and the results have been averaged over 2 • 10 4 and 4 • 10 3 trajectories respectively. For a given value of p, 
we have performed a series of simulations for different values of the mass of the tagged particle. Taking due account 
of the constraint cq < 1, the value of the tagged particle's mass must be smaller than 10 2 for p = 0.1 and 10 3 for 
p = 0.01. 

Figure [T] shows the evolution of the ratio T/T g for a system with p = 0.01 and for two values of the tagged particle's 
mass, that is, m — 20m g and m = 200m 9 , as a function of the number of collisions per particle, t, defined as 



r = 



dt'Ug(t') 



9 j.* 



2(1 - Co) le 



(61) 



The initial value of temperature of the tagged particle is T(0) = for all the trajectories, since at t = 0, the intruder 
has a prescribed velocity. As we can see, in this scale, the evolution to the stationary value of the ratio of the 
temperatures is slower as we increase the mass of the tagged particle. This implies that there are values of the tagged 
particle's mass for which the ratio of the temperatures will not reach its stationary value in the time of the simulation 
(for instance, the number of particles in the bath for r = 200 is N g ~ 1800, which hinders correct statistical sampling). 
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Figure 4: (Color online) Reduced fourth moment of the tagged particle velocity as a function of the dimensionless time r for 
p — 0.1 and different values of the t 



Figure 5: (Color online) Mean velocity of the tagged particle as a function of t*/(l — eo) for a system with p = 0.1 and several 
values of the tagged particle's mass. The dashed line is the theoretical prediction given by Eq. 1|33[) . Inset: Mean velocity of 
the tagged particle as a function of t* /(l — eo) for a system with p = 0.1 on a logarithmic scale. 

The theoretical prediction in this time scale (dashed line) is obtained directly from Eq. (|34|>. taking into account Eq. 
(|6ip . The agreement between theory and simulations is good. 



Similar simulations were performed with different values of the tagged particle mass (m ranging from 15 to 100 for 
p = 0.1 and from 15 to 900 for p = 10" 2 ). Since T g (Q) = 0, Eq. ((34]) predicts 



As can be observed in Fig. [21 all simulation data for $ collapse onto a single curve. The time scale t* can be calculated 
from the scale r defined in (|61|) . In the same vein, we can obtain the stationary value of the temperature ratio, which 
is plotted as a function of A = m g /m in Fig. [3]; this validates our theoretical analysis. 

In order to probe -at least partially- the Gaussian nature of the time dependent tagged particle velocity statistics, 
we have measured the reduced fourth moment 4(w 4 ) / (d(d + 2)(v 2 } 2 ). As can be seen in Fig. [5] where we have plotted 
the results for p = 0.1 and several values of A as a function of the number of collisions per particle, r, the value of 
this quantity is in agreement with the Gaussian prediction (that is unity) within the statistical uncertainties. 

Consider next the mean velocity of the tagged particle, u(t*). In order to study the decay of this quantity, we have 
performed a set of simulations starting with a component of the velocity field in the x direction when it is immersed 
in a bath in the homogeneous decay state, u x (0). In Fig. [5J we plot u x (t*)/u x (0) as a function of t*/(l — eo) for 
p = 0.1 and several values of m. In this time scale, the data for all values of m collapse due to Eq. (f33|) . In the inset, 
the same quantity is plotted on a logarithmic scale. If the theoretical prediction in Eq. (|33| is verified, the above 
plot must lead to a straight line with slope £ u = 1 (dashed line), where Cu is the decay rate of the mean velocity. On 
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Figure 6: Decay rate of the mean velocity as a function of A for a system with p = 0.1. The symbols are from DSMC simulations 
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dashed line is the theoretical prediction, Eq. (I53|l . Inset : same quantity as a function of n g (0) 2 /n g (t*) 2 . 
the theoretical prediction of Eq. (|63|l where B follows from (|65|l . 



The 

The dashed line is 



the other hand, £ u can be fitted on the logarithmic plot of the inset. Reporting the corresponding measures in Fig. 
[5] against the mass ratio, it appears that the theoretical prediction £ u = 1 is approached as A — > 0, as expected. 

Finally, the accuracy of the prediction for the diffusion equation has also been tested by measuring in DSMC 
simulations the mean square displacement of the tagged particle. In contrast to the granular case phenomenology 
and as a consequence of the continuous decay of particles in the bath, the mean squared displacement increases 
exponentially in the t* scale, see Eq. (I5B")) . In Fig. [7] we have plotted the time evolution of (r 2 ) in the scale t* for 
a system with p = 0.01 and m = 60m g . The dashed line is the theoretical prediction given by Eq. (|53p . and shows 
good agreement with numerical data. The same quantity, written in terms of the bath density, reads 



<r 2 (0> = dD £ 2 (0)- 



D 



nM 2 

n g (t*y 



- 1 



(63) 



where we have defined 



B = 



dD £ 2 {0) 



(64) 



and we have taken into account that n g (t* 



,(0)e 



-e*i* 



It then appears that the mean squared displacement 



increases linearly with n g (0) 2 /n g (t*) 2 . This is full agreement with the simulation results, see the inset of FigEl Such 
a plot allows us to extract by linear fitting the coefficient B, which can then be compared against the prediction of 
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t igure 8: Values of Bn 2 g (0) as a function of A for systems with p = 0.1 (left) and p — 0.01 (right). The dashed line is for the 
prediction of Eq. (|65|l . 



Eq. which explicitly reads 



16 A /2d 2 7r 1 - d ^ (1 - d) r (I) 



B = 



n2(0)p(-16 + V2£)(16- 



y/2(4d-l)%) 



(65) 



where we have taken into account that uq — a g . Such a comparison is worked out in Figure [8] and fully corroborates 
the theoretical analysis, with again an improved agreement when A decreases. 

It would be also interesting to confirm our theoretical predictions with Molecular Dynamics simulations. In the low 
density limit, it is expected to get similar results. In fact, some simulations were performed finding qualitatively the 
same behavior but with much more statistical inaccuracies. 



V. CONCLUSIONS 



In this paper, the diffusive behavior of a tagged intruder immersed in a gas of particles undergoing ballistic anni- 
hilation (i.e. which annihilate with probability p or scatter elastically otherwise), has been analyzed. The collisions 
between the tagged particle and the surrounding gas are elastic. Some similarities are found between our system 
and the elastic or inelastic case [HI, [l8|], but, on the other hand, important differences arise as a consequence of the 
continuous decay of particle number in the system. 

We start from the Boltzmann-Lorentz equation for the distribution function of the tagged particle, which is valid, 
in principle, for arbitrary mass of the tagged particle. In the limit of a very massive tagged particle, a Fokker-Planck 
equation for the distribution function is derived by means of a systematic expansion in the mass ratio A. Our approach 
holds in the limit A « 1, but we additionally have the more stringent condition that the parameter introduced in 
Eq. (|2"5|) , e cx j^, must be smaller than unity. Analysis of the Fokker-Planck equation leads to predictions for the 
temperature ratio, the decay rate of the mean velocity of the tagged particle and the diffusion coefficient. As in the 
inelastic case 18j, the theory predicts that the ratio between the temperatures of the tagged particle and the gas is 
constant in the long time limit as a consequence of equilibrating cooling rates. When represented in the appropriate 
time scale, which is proportional to the number of collisions experienced by the tagged particle, temperature ratios 
collapse for all values of A and p considered. Likewise, the mean intruder velocity (averaged over bath realizations) 
decays exponentially. The dynamics of the distribution function of the tagged particle is governed by a Fokker-Planck 
operator which spectral properties are known. More specifically, as the eigenvalues of this operator are non-positive, 
the distribution of the tagged particle approaches a Gaussian in the long-time limit and admits a scaling form similar 
to the one for the distribution function for the particles in the gas but with a different temperature. At variance with 
the situation of an elastic intruder in a bath of inelastic particles [22| , there is apparently no mapping between our 
problem and a well chosen elastic system. A unique vanishing eigenvalue is responsible for the slow diffusive behavior 
of the tagged particle density. The corresponding diffusion equation has been derived in the hydrodynamic limit, by 
means of a projector decomposition, which yields an explicit expression for the diffusion coefficient. From a different 
point of view, the expression for the mean squared displacement has also been derived "a la Einstein" . Following 
this route, the diffusion coefficient is expressed as a Green-Kubo formula in terms of a weighted time integral of the 
tagged particle velocity correlation function. This provides a more physical perspective on the results derived from the 
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projector method. As already mentioned, the mean squared displacement for this system does not increase linearly 
in the collision per particle time scale, as is the case in the elastic and inelastic cases. This different behavior is due 
to the time dependent bath density. In the elastic case, both the temperature and the density do not depend on 
time. In an inelastic system [3, the time dependence goes through the temperature and could be absorbed in 
the collision per particle time scale, which turns out to be impossible in our system where the mean free path £(t*) is 
an increasing function of time. 

Finally, our analytical results have been tested by numerical simulations, and a very good agreement has been 
reported for all the range of parameters considered. As expected, the agreement is all the better as A is smaller. 
In summary, the work reported here provides an example of the accuracy of hydrodynamics to describe a system in 
which there are no conserved quantities in binary encounters (no collisional invariants). 
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Appendix A: SOME USEFUL APPROXIMATIONS 



For the sake of completeness, we provide here the approximate expressions for the density and temperature decay 
rates, that are relevant for explicit computation of several of the quantities discussed in the main text. They have 
been obtained from a truncated Sonine expansion (Sonine polynomials being particular types of Laguerre polynomials, 
particularly convenient for kinetic theory calculus) 0, [g] . 

C, = ^fl-a 2 AV (AD 



where 



4 V 16 , 

d + 2 / 8d+ll\ 
Ct = ^ l + a 2 ^-), (A2) 



8(3-2V2> 

CL2 = 7= 7= ■ (A3) 

{4d + 6 - V2)p + 8V2(d - - p) 



Appendix B: FROM THE BOLTZMANN-LORENTZ EQUATION TO THE FOKKER-PLANCK 

EQUATION 

In this Appendix we expand the collision operator, Eq. @, in series of A. We start by multiplying the collision 
operator by a generic function H(v) and integrate in velocity space 

fdvH(v)J[r,v,t\FJ] 

= at 1 |dv|dv 1 i/(v)|d«Te(g ■ <x)(g • &)[F(v*)f(vl) - F(v)/( Vl )]. 



The above expression can be written 

fdvH(v)J[r,v,t\F,f] 



at 1 Jdyfdw 1 F{w)f{v 1 ) Jd&Q(g ■ <x)(g • &)[H(y - <5v) - H(v)}, 



(Bl) 



(B2) 
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where we have introduced 



r 2A 

dv = —(& ■ cr)(J, 

1 + A vs ; , 



(B3) 



which is the increment of the tagged particle velocity due to collisions with a particle of the bath (it should be 
remembered that g = v — vi). Equation (|B2[) essentially tells us how the function H varies due to collisions. If we 
admit that A is small enough, we can expand H{w — <5v) around v in powers of 5v, keeping only the lower orders 



dH(v) 
<9v 



H (v - 5v) ~ H(v) - 
If we introduce expansion (|B4[) in equation (|B2|) we obtain 

fdvH(v)J[r,v,t\F,f] 



1 



where we have introduced 



dv#(v) ^ — • [A(v)F(v)] 



A(v) = 



(i + A)r(^) 



ld_d_ 
2 <9v dv 



[N(v)F(y)] 



(B4) 



(B5) 



(B6) 



N(y) = 



2A 



1 



7T 2 (T 



d-1 



dv i/( v i) 



cZ + 3 o T 3 

^T 5l+ 2- 9 



3 3 I 



(B7) 



In the last expression I is the unit tensor. As H (v) is a generic function of v, we can compare equations (|B1[) and 
(|B5|) . and we obtain that the collision operator can be written as 



J[r,v,t\FJ] ~ A . [A(v)F(v)] + ±AA : [JV(v)F(v)]. 



(B8) 



We next specify A and N within the scaling form provided by the homogeneous decay state of the bath. This will 
lead us to identify the remaining A dependence in these coefficients and to simplify the functional dependence in the 
tagged particle velocity. To this end, we introduce the dimensionless velocities 



/ 2T (t)\^^ /2T(t)\ L ^ z 

where v g (t) = ( ^ ' \ and v (t) = I ^' \ , with T(t) the temperature of the tagged particle and T g (t) the 
temperature of the suspending gas. The relative velocity g = v — vi can be written as 



vo(t)' 

1/2 



Cl 



%(*)' 



(B9) 



g = Vg(t) 



T(t) 
T g (t) 



1/2 



A 1/2 c- u s (t)ci. 



A formal expansion in (T/T g )A leads to 

A(v,t)= 7 (i)v, N(y,t)=2j(t)I. 
The definitions of 7 and 7 are respectively 

7 (i) = j e [n g (t),T g (t)]a(p), 



l(t) = >y e [n g (t),T g (t)]a(p)b(p) 



Tg{t) 



(BIO) 

(Bll) 

(B12) 
(B13) 



where ~/ e (t) is the same friction coefficient as for an elastic system at the corresponding density and temperature 



le[n g (t),T g (t)] = 



4-7T ~ 

dT(d/2) 



& 1/2 n g (t) 



2T g (t)\ 1/2 , 1 



rn 



'0 > 



(B14) 
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and a, b are functionals of the distribution function of the bath which depend only on the parameter p 



a(p) 



T(d/2) 



r((d- 



2) f 

jyTjj Jdc lX H{ci)c u (B15) 
= * I^MfOff , (B16) 
These coefficients have been evaluated in the first Sonine approximation and depend very weakly on p 

a( P ) = ^f^, (B17) 

= 8 + 3^ 
8 - a 2 (p) 

where a 2 , defined in (|A3|) . is the gas velocity distribution kurtosis (a Gaussian ansatz would amount to setting a 2 = 0). 
By dimensional analysis and taking into account the explicit formulas for A and N, equation (|B11[) . we can see that 
the two terms we have considered in the expansion of the collision operator, Eq. (|B8[) . are of order n g v g a ~ l A, while 
the other terms in the Kramers-Moyal expansion are at least of order n g v g a'^~ 1 A 3 / 2 . Hence, we can conclude that 
the leading order contribution in A of the collision operator is actually the one written in (|B8[) . 



Appendix C: EQUATIONS FOR THE VELOCITY AND TEMPERATURE OF THE TAGGED PARTICLE 

In this Appendix we derive the equations for the mean velocity and temperature of the tagged particle. Taking 
moments in the Fokker-Planck equation, Eq. (fT5| . we obtain for the velocity 



du(t) 



()t 



a(p)b(p)% e (t)-^\ F(r,v,t). (CD 



By integration we have 



Jdrjdvv^-(v V) +7e(*)a(p)^ ■ v 

T d 2 
m ov z 

^ = JdrJdvv-y e (t)a(p)^.(vF(r,v,t)) (C2) 

= - JdrJdvi e {t)a(p)vF(r,v,t) (C3) 

= -7e(*)a(p)u(i). (C4) 
Consequently, the equation for the mean velocity is 

du(t) 
dt 

Taking into account the definition of temperature, 



- le (t)a(p)u(t). (C5) 



~T(t) = Jdrjdv^m[w^u(t)] 2 F(r^,t) 



drldv-m[v 2 - u 2 (t)]F(r,v,t). (C6) 



2 

we can write 

d f, o^/ \ . . 9u(t) 



2 dt w 2 



dvv 2 F(r, v, i) - 2u(f) 



dtj j v ' w w at 

In order to evaluate the first term on the right hand side, we make use of the Fokker-Planck equation: 



(C7) 



-| jdvfdwv 2 F{v,w,t) = -2 7e (t)a(p) [drfdvv 2 F{r,v,t) + ^T gle {t)a{p)b{p). (C8) 

r elocity into account, we obtain 

= -d le (t)a(p)[T(t)-b(p)T g (t)}. (C9) 



Taking this formula and the equation for the velocity into account, we obtain 

ddT(t) 



2 dt 
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Appendix D: THE DIFFUSION EQUATION 



In this Appendix we derive the diffusion equation for the tagged particle's density. The starting point is the 
Fokker-Planck equation (I43[) 



Ap k (v*,t*) = [A FP (v*) - i£ (t*)k • v*]F k (v*,**), 



in which we introduce the two projectors 



P9(v*) = <XAf(v*)| 3 (v*)) XM (v*), 
P ± g(v*) = (l-P)s(v*). 



(Dl) 



(D2) 
(D3) 



Here, we have introduced the maxwellian distribution, xm(v*), which is the eigenfunction of App associated with the 
eigenvalue and we have used the scalar product defined as 

(/(v*)| 5 (v*)) = J <fr*XMV)/V)0(v*), (D4) 

f> being the complex conjugate of /. In a next step, we decompose the function P k in PP k and Pj_P k , and write the 
equations for these two quantities 



JL+U Q (t*)P\L-v* -PA FP 



PF k = -^ (t*)Pk-v*PxF k , 



(D5) 



dt* 



i£ (t*)P±k-v* -P ± A FP 



P±F k = -i£ (t*)P 1 _-k-v*PF k . 



(D6) 



We are interested in obtaining a closed equation for PF k in the hydrodynamic limit. To achieve this goal, we formally 
solve the equation for Pj_F k 



P±P k (vV*) = G (t*)F k (v*,0)- f dt*'G t .,(t*-t*')P ± i£ (t*')k-v*PF k (v*,t*'), 

Jo 



where the operator Gt*'(t* — t*') is defined as 

d 



G t „(t* - t*') = P ± [A FP (v*) - i£ (t*)k ■ v*]P ± G t „{t* - t*'), 



with G t »(0) = 1. In the long time limit, the term associated to the initial condition vanishes and we have 

P±P k (v*,t*) = - / dt*G t *>(t* -t*')P x i£ (t* - £*')k-v* PP k (v*, t*-t*'). 
Jo 



(D7) 



(D8) 



(D9) 



In order to obtain the diffusion equation to order fc 2 , we only need Pj_P k to order fc, so we write the expression for 
Gf{t* — t*') to leading order 



G t .- t „(t*')~e 



P ± A FP P ± t" 



(D10) 



We then have 



Pj_F k (v*,f) ~ -/ dt*'e p±AFpP±t " P ± i£ (t* ~t*')k- v*PF k (v*,i* -t*') 
Jo 

~ -£ (t*) f dt*'e p±AFpP±t "- t " t " P^ik- v*PF k (v*,i* -t*'), 



(Dll) 



where we have used that £o{t*) ~ e e ***. We subsequently have to relate PP k (v*,i* — t*') with PF k (v*,f )■ To be 
consistent with the hydrodynamic approximation, this is done to leading order 



PF k (v*,t* -t*')~e 



-PA FP Pt" 



PP k (v*,i*) = (XM(v*)|F k (v*,i*)) XM (v*). 



(D12) 



17 

Now we can write the equation for the Fourier transform of the density, rik(i*) 

«k(**) = -lke (t*)} 2 ^ J dv*v* • J rfr / e p - A - p - t *'e- e ***'P ± v* X M(v*)n k (i*). (D13) 



<9i 

In other words, 

^«k(f) = -D [£ (t*)k] 2 n k (t*). (D14) 

where 

D Q = i y dv*v* • jf dr'e( Pj - A ^ p ^- £ *) t *'P J _v*XM(v*). (D15) 
Finally, we can evaluate Do exactly since v|xm(v*) is an eigenfunction of Kpp with eigenvalue Ai = — 1 

D = J dv*v* J* dt'e^-^'P^XM^*) = (D16) 
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